The Einstein static universe in Loop Quantum 
Cosmology 



Luca Parisi 1 ' 2,3 , Marco Bruni 3,4 ' 5 , Roy Maartens 3 
and Kevin Vandersloot 3 

1 Universita di Salerno, 84081 Baronissi (Salerno), Italy 
2 INFN, Sezione di Napoli, Salerno, Italy 

institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth POl 
2EG, UK 

4 Dipartimento di Fisica, Universita di Roma "Tor Vergata", 00133 Roma, Italy 
5 INFN, Sezione di Roma "Tor Vergata", 00133 Roma, Italy 

Abstract. 

Loop Quantum Cosmology strongly modifies the high-energy dynamics of 
Friedman-Robertson- Walker models and removes the big-bang singularity. We 
investigate how LQC corrections affect the stability properties of the Einstein static 
universe. In General Relativity, the Einstein static model with positive cosmological 
constant A is unstable to homogeneous perturbations. Using dynamical systems 
methods, we show that LQC modifications can lead to an Einstein static model which 
is neutrally stable for a large enough positive value of A. 
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1. Introduction 

The Einstein static universe in General Relativity (GR) is a closed Friedman- Robert son- 
Walker (FRW) model that is unstable to homogeneous perturbations pQ. (Note that it 
is neutrally stable to inhomogeneous scalar perturbations with high enough sound speed 
and to vector and tensor perturbations [2].) 

The stability of Einstein static models in high-energy modifications of GR is an 
interesting mathematical question, and is also relevant for scenarios in which the Einstein 
static is an initial state for a past-eternal inflationary cosmology, the so-called Emergent 
Universe scenario [3]. The standard model of cosmology (see e.g. [I]) has a flat infinite 
spatial geometry, and experiences a period of primordial inflationary expansion, which 
is preceded by a big bang singularity in the classical theory. Observations however do 
not prove that the geometry is flat: the Universe could have nonzero spatial curvature, 
as long as the late-time effect of this curvature is very small. In particular, a positive 
curvature allows for an "Emergent Universe" that originates asymptotically in the past 
as an Einstein static universe, and then inflates and later reheats to a hot big bang era. 
This model generalizes the Eddington-Lemaitre model [5], and is a counter-example 
to the notion that inflation can never be past-eternal and thus cannot avoid an initial 
singularity - because it is closed, it avoids the theorems showing that inflation cannot 
be past eternal [6]. 



Generalizations of the Einstein static solution in high-energy modifications to GR 
have been considered in the Randall-Sundrum braneworld scenario [7j and in f(R) 
theories [8]. The homogeneous stability of Einstein static models in R + aR 2 gravity 
has also been investigated [9J. 

Another theory leading to high-energy modifications of GR is Loop Quantum 
Cosmology (LQC) [10] , which is a canonical quantization of homogeneous cosmological 
spacetimes based on Loop Quantum Gravity [TT] . The gravitational phase-space 
variables are an su(2) valued connection and conjugate triad, and the elementary 
variables underlying the quantization are the holonomies of the connection and the fluxes 
of the triad. The quantum theory obtained from LQC turns out to be inequivalent 
to Wheeler-de Witt quantization (the LQC polymer representation is different from 
the usual Wheeler-de Witt Schrodinger representation). Wheeler-de Witt quantization 
does not resolve the cosmological singularity, but in LQC a generic resolution of such 
singularities has been obtained. Initial quantizations of LQC lead to a regularization 
of the big bang singularity [12] resulting from the fact that the quantum Einstein 
equation is non-singular as well as from modifications to the scalar field energy density 
and dynamics. The modifications to the scalar field dynamics were based on effects 
arising from quantum inverse scale factor operators. Subsequently, the elucidation of 
the consequences of using holonomies as the basic variables has shown that gravity 
is modified [13], i.e., the structure of the Friedman equation, rather than solely the 
energy density of the matter field. These gravitational modifications typically become 
important at lower energy scales than the modifications to the scalar field dynamics. 
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Mulryne et al. [2] used the scalar-field modification approach to investigate the 
stability of the Einstein static model to homogeneous perturbations. They found that 
the new LQC Einstein static model is a centre fixed point in phase space, i.e. a neutrally 
stable point, for a massless scalar field with w = p^/p^ = 1. This modification of 
stability behaviour has important consequences for the Emergent Universe scenario, 
since it ameliorates the fine-tuning that arises from the fact that the Einstein static is 
an unstable saddle in GR. 

Here we consider the same question, but using the LQC gravitational modifications, 
and neglecting the higher energy modifications to matter. We consider a perfect fluid 
with p = wp and w > — 1. We show that, for all w > —1/3, there is a new centre, i.e. a 
neutrally stable point representing an Einstein static model, but only if the cosmological 
constant A is above a critical scale, 

A>6.6vrMj,. (1) 
2. Critical points of closed FRW models in LQC 

The loop quantum effects that we investigate manifest themselves in the form of a 
modification to the classical Friedmann equation. For the closed FRW model, the 
explicit form is given by [13] 

3 3 a z J \ p c np c Kp c a\' 

where n = 8nG = 8n/Mp, and the critical LQC energy density is 

p c « 0.82M£ . (3) 

It is evident that the first term in parentheses is the classical right hand side of the 
Friedmann equation, with the quantum modifications appearing in the second term. 
The classical GR limit is achieved in the limit as p c goes to infinity whence the second 
term approaches unity. The classical energy conservation equation continues to holcj||, 

p = -3Hp (1 + w) . (4) 

Note that if 2 > imposes the limits 
3 3 

— < up + A < k Pc + — . (5) 

cr er 

The modified Raychaudhuri equation follows from Eqs. (@J and (j2J): 

K , , / 2p 2A\ 
H = --p (1 + w) ( 1 - — - — ) + 

— 6 ( 6 ) 

a? Kp c a A 



2 V Pc K PcJ 

+ 



^ 2p 2 A 3p(l + w) 



Pc Kp c p c 



\ Note that this means we are not considering the inverse volume effects that would modify the scalar 
field energy density and hence modify the energy conservation equation. For the closed model there 
are indications that the inverse volume effects are negligible if it is required that the universe reach 
macroscopic size [13) . 
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We will find the critical points to the system of Eqs. (jlj), (jB]) as well as 

a = aH , (7) 

which follows from the definition of H. The solution space is a 2-dimensional surface 

in the three-dimensional (p,a,H) space, defined by the Friedman constraint (J2J). The 

system admits two critical points, which are static solutions a = H = p = 0. The first 

critical point is the standard GR Einstein static universe, while the second is a new 

LQC Einstein static universe: 

2A 2 2 
Pgr - n , o \ > a GR ~ J^~, s > ( 8 ) 

_ 2 (A - k Pc ) 2 2 

~~ «(1 + 3w) ' a ^- K p iQ (l + W )- (J) 
The conditions under which these static solutions exist are summarized in Table [U and 
follow from a 2 , p > 0. 

A remarkable feature of the new LQ fixed point is that it is possible to have an 
Einstein static universe even for vanishing cosmological constant. Indeed, as one can 
see from Eq. (Q and Table [21 when A = the LQ fixed point exists and is unstable. 

Table 1. Conditions for the existence of Einstein static critical points. 





GR 


A > 


w > -1/3 






A < 


-1 < w < 


-1/3 


LQ 


A < Kp c 


-1 < w < 


-1/3 




A > Kp c 


w > -1/3 





For the system of equations (BJ, (jSJ) and (JTj), linearized stability analysis fails to 
give complete information about the properties of the two critical points, when they 
both exist, since in all cases, for one of the two points there is always a pair of complex 
eigenvalues with vanishing real part. In this case the linearization theorem (see e.g. [15] ) 
tells us that a fixed point which is a centre for the linearized system is not necessarily a 
centre for the full nonlinear system. In addition, the linearization of equations (HI), (EI) 
and (j7]) always leads to a third vanishing eigenvalue, simply because the actual dynamics 
is two-dimensional because of the modified Friedman equation (j2J). Therefore, using this 
constraint, it is convenient to rewrite the system reducing the number of equations to 
two, the dimension of the configuration space. Solving Eq. (j2J) for a 2 , we find that 

a 2 = f ± (p,H), (10) 



where 



3 



2{Kp + A) + k Pc (l ± y/l - \2H 2 /k Pc 



f±=«—, . 7T2 , ' Tn" ■ (11) 



2 ( K p + A)' + k Pc (3H 2 - up- A) 
Substituting this into Eq. we find two branches for the time derivative of the Hubble 
parameter, 

LQ : p = -3Hp (1 + w) and H = F + (p,H), (12) 
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GR:p = -3Hp{l + w) and H = F„{p,H), (13) 

where 



k,, , A 2p 2A . 



2 V Pc ^Pc/ «Pc/± 

1 

- 1 - 

Pc UPc ' Pc_ 



f: 



± 



(14) 



In the classical limit, 

lim/ + = 0, (15) 

lim /- = ^m~J ' ( 16 ) 

where the second equation is the GR Friedman equation. This shows how the two 
branches in Eq. ffTUl) recover the GR and the new quantum static solution. In addition, 

lim F + = oo, (17) 

lim F_ = -H 2 - ^(1 + w) + ^ , (18) 

where the first limit is consistent with Eq. (1151) . and the second limit gives the GR 
Raychaudhuri equation. 

The system ( 1T2"]) admits the static solution with p L Q as in Eq. Substituting this 
into Eq. (fTTj) reproduces o\q as in Eq. dHJ). We evaluate the eigenvalues of the Jacobian 
matrix at this point, to find the two eigenvalues 

\ LQ = ± v /(np c -A)(l + w). (19) 

Thus the LQ fixed point is either unstable (of the saddle kind), when np c > A and 
— 1 < w < —1/3, or a centre for the linearized system, i.e. a neutrally stable fixed 
point, when Kp c < A and w > —1/3. (The limits take into account the conditions in 
Table [TJ) Again, for the latter point the linearized analysis is not sufficient and therefore 
we turn to a numerical analysis in the next section. 

For the system Eq. fflBl . we find the GR static solution, and the eigenvalues of the 
linearized system are 



\ GR = ±^A(l + w). (20) 

These are real with opposite signs for A > and w > —1/3, so that the fixed point 
is unstable (of the saddle type). For A < and —1 < w < —1/3, the fixed point is a 
centre, as confirmed by numerical analysis. It does not appear to be widely known that 
the GR Einstein static universe can be neutrally stable when A < (see also [9]). 
The results of the linearized stability analysis are summarized in Table |2j 



3. Numerical integration 



In order to extend the linearized stability analysis, we perform numerical integrations 
of the systems ( TT2"j) and (ITBl . We also integrate the nonlinear system of Eqs. 01]), fl§D 
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Table 2. Eigenvalues for the critical points in Table Q] 









Ai 




A 2 


GR 












A > 


and 


w > -1/3 


> 




< 


A < 


and 


-1< w < -1/3 


Re(Ai) 


= 


Re(A 2 ) = 








Im(Ai) 


> 


Im(A 2 ) < 


LQ 












A < np c 


and 


-1 < w < -1/3 


> 




< 


A > np c 


and 


w > -1/3 


Re(Ai) 


= 


Re(A 2 ) = 








Im(Ai) 


> 


Im(A 2 ) < 



and (J7J), with initial conditions fulfilling the Friedman constraint (|2J), in order to show 
the full configuration space diagrams and gain a better understanding of the dynamics. 
We take A > and w — 1. 

First we consider the case A > np c . The (H,p) plots in Figs. [Hand [2] are obtained 
by integrating the system f[T2"j) and (TX5T) in some neighbourhood of the fixed points. The 
plots show that the GR Einstein static solution is a saddle while the LQ Einstein static 
solution is a centre. 




Figure 1. Dynamical behaviour of the system around the GR fixed point for the case 
A > np c , with A/k = 2, w = 1 (using units Mp = 1). 

A better understanding can be obtained by plotting the whole 3D space (H, p, a) 
for a wide range of initial conditions, shown in Fig. The trajectories lie on the 
2-dimensional Friedman constraint surface. 

The behaviour near the fixed point is in agreement with the linearized stability 
analysis, but new interesting features arise. For initial conditions far enough from the 
fixed point, there are trajectories that wrap around the Friedman tube, so that cyclic 
models are possible even if they are not related with the centre fixed point, since they 




Figure 2. As in Fig. [I] for the LQ fixed point. 




Figure 3. Trajectories on the Friedman constraint surface, for the same parameters 
as in Figs. [T] and [2] The GR fixed point is at the bottom of the Friedman surface, 
while the LQ point is at the top. Note that some trajectories wrap around the "tube" 
but cannot be shrunk continuously to the LQ fixed point. 



cannot be shrunk to a point; see Fig. HI 

This behaviour can also be interpreted via a plot of the energy density against a. 
Defining 

3 3 A 

P- = o j P+ = Pc H o , Pm= P + -, (21) 



Eq. ([2]) can be written in the form 



H 2 = ^~ ( Pm ~ P-) {P+ ~ Pm) ■ (22) 
dpc 

When p m = p_(a) or p+(a), then H—0, so the system undergoes a bounce or starts a 
recollapsing phase, respectively. This is illustrated in Fig. [5j 
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Figure 4. As in Fig.[3l showing trajectories that start far from the fixed points. 



p 




0.16 0.18 0.2 0.22 0.2A 0.20 0.20 0.3 0.32 



Figure 5. The upper red (grey) curve is p+(a) and p-(a) is the lower curve. 
Trajectories for several initial conditions are depicted in black. The GR fixed point is 
the box on the curve p_, and the LQ point is the circle on the curve p+. 



The second case A < np c is illustrated in Fig. [61 For w — 1, only the GR fixed point 
is present and it is unstable. The trajectories are again wrapped around the Friedman 
surface. During a contracting phase they tend to a minimum of the scale factor and 
then, after a bounce, there is a phase of expansion. 

We can integrate Eq. (0J for p as a function of a, and substitute the expression for 
a(p) into the modified Friedman equation ([2]). Using the dimensionless variables, 



2 3k/9 c 2 up np c 

x =^ H ' y = T' 5 = X 



this leads to 



x 2 



\j + 1 - Cy 2 ^ w+ ^] [B-y-l + Cy 2 ^ w+ V] , (23) 
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Figure 6. The case A/n = 4 x 10 7 , with w = 1. 

where C is a constant of integration. This expression allows us to get a better 
understanding of the particular case depicted in Figs. [3] and HI where both the fixed 
points are present, by finding the separatrix. The separatrix is a curve (actually, a 
union of orbits) which marks the boundary of regions where the dynamical behaviour 
of the system is different. In this case the separatrix is the junction of the stable and 
unstable manifold with the GR Einstein static hyperbolic fixed point. First we solve 
Eq. (1231) for C, then we substitute numerical values for the parameters w = 1 and 
B = .4121769562 (using M P = 1), and we also substitute x = and y = 1/2 for 
the GR fixed point. This produces a numerical value for C . Through this procedure, 
the equation of the separatrix is implicitly given by Eq. (123]) with fixed values of the 
parameters C, w and B. 

More insight can be acquired via the plot of the separatrix projected onto the (x, y)- 
plane, as shown in Fig. [7J The curves depicted are the separatrix, which wraps around 
the Friedman surface, and some other trajectories which are first integrals obtained for 
different values of the integration constants. The closed loops around the Friedman 
surface are outside the region marked by the separatrix; the other curve is not closed 
around the tube, and it can be shrunk continuously to the LQ fixed point. 

4. Conclusions 

We have shown that LQC gravitational modifications to the Friedman equation lead to 
a new high-energy critical point for the Einstein static universe. This LQC Einstein 
static model is an unstable saddle (like the standard GR Einstein static solution) for 
sufficiently negative pressure and a sub-critical cosmological constant, i.e., 

-l<w<--, A<k Pc . (24) 
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Figure 7. Projection of the phase space onto the (x, y)-plane for the case of Figs. [3] 
and [4j It represents the trajectories as seen from the bottom of the Friedman surface. 
Both the GR fixed point (a box) and the LQ fixed point (a circle) exist. The two 
vertical lines are the edges of the Friedman surface. The red curve represents the two 
branches of the separatrix. 



If w is large enough and the cosmological constant is above the critical value, then the 
LQC Einstein static model is a centre: 

w>--, A>np c . (25) 

This neutrally stable behaviour is in strong contrast to the GR case, where the Einstein 
static model is unstable for all (positive) values of A. 

Modified stability of the Einstein static model is also found for R + aR 2 
gravity [9j and for the LQC case with matter modifications, rather than gravitational 
modifications [14] . This illustrates the general point that high-energy modifications to 
GR, which typically strongly modify the big bang singularity, also modify the dynamical 
nature of the non-singular Einstein static universe. 

Our result means that the fine-tuning problem for the Emergent Universe scenario 
is qualitatively changed by LQC gravitational modifications. In GR, the Einstein static 
is an unstable saddle, so that severe fine-tuning is required if the Einstein static is to 
be the initial state for a past-eternal inflationary cosmology. With LQC gravitational 
modifications, the initial state becomes a centre when Eq. ( |25i) holds, thus ameliorating 
the fine-tuning. The same point was made in Mulryne et al. [H] in the case of LQC 
matter modifications. In that case however, the analysis is restricted to a massless 
scalar field (w = 1), and the centre exists for all A. Our analysis applies for all perfect 
fluids with w > — 1, since we focus on the LQC gravitational modifications. These 
gravitational modifications operate typically at a lower energy scale than the matter 
modifications. Clearly, from the standpoint of the Emergent Universe scenario, one 
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needs a further mechanism to break the infinite cycles about the centre fixed point, 
leading to a subsequent inflationary phase as discussed in [16J for a general class of 
braneworld models sourced by a scalar field with a constant potential. This mechanism 
will be the subject of future work. 
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